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ABSTRACT 

We report the first astrometrically determined mass of an extrasolar planet, 
a companion previously detected by Doppler spectroscopy. Radial velocities first 
provided an ephemeris with which to schedule a significant fraction of the HST 
observations near companion peri- and apastron. The astrometry residuals at 
these orbital phases exhibit a systematic deviation consistent with a perturba- 
tion due to a planetary mass companion. Combining HST astrometry with radial 
velocities, we solve for the proper motion, parallax, perturbation size, inclination, 
and position angle of the line of nodes, while constraining period, velocity ampli- 
tude, longitude of periastron, and eccentricity to values determined from radial 
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velocities. We find a perturbation semi-major axis and inclination, a = 0.25 ± 
0.06 mas, i = 84°±6°, and Gl 876 absolute parallax, n abs = 214. 6± 0.2 mas. 
Assuming that the mass of the primary star is = O.32M , we find the mass 
of the planet, Gl 876b, M b = 1.89 ± 0.34M Jnp . 

Subject headings: astrometry, stars: distances, planetary systems, techniques: 
interferometric 



1. Introduction 

In 1998 two groups searching for unseen companions to stars using the technique of 
Doppler spectroscopy announced the discovery of a companion to the M4 dwarf star, Ross 
780 = Gl 876 (Delfosse et al. 1998, Marcy et al. 1998). This companion was characterized by 
Msini ~ 2M Jup and P~60 d . Later, a second companion was detected, Gl 876c (Marcy et al. 
2001), with Msini = 0.5QM Jup and P~30 d . We were granted HST time to characterize the 
perturbation due to the longer period companion, GL 876b. The low mass of the primary 
(M* = 0.32M Q ), the period of Gl 876b, and the system proximity (4.7pc) suggested that 
even an edge-on orientation (M = Msini) would produce a detectable perturbation, thus 
a companion mass. A transiting companion of HD 209458 allowed Henry et al. (2000) to 
establish the first mass of an extrasolar planet. We report the first astrometrically determined 
mass of an extrasolar planet. 



2. Observations and Reductions 

Observations were obtained with FGS 3 in fringe tracking mode. The FGS instrument 
is described in Nelan & Makidon (2002). Each of our 27 data sets required ~ 40 minutes of 
spacecraft time. Data were reduced and calibrated as detailed in Benedict et al. (2002b), 
McArthur et al. (2001), and Benedict et al. (1999). While acquiring each set, we measured 
reference stars and the Gl 876 multiple times to correct for intra-orbit drift (e.g. Benedict 
et al. 2002a, Fig. 1). Three prime epochs (1, 3, and 4, boldface in Table 1) contain 30 or 45 
observations of Gl 876 and were obtained near companion peri- and apastron, as predicted 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555 
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from the radial velocity (RV) data. For a planetary mass companion of unknown inclination 
observations with HST/FGS precision at any other orbital phase are less likely to yield a 
detection. To minimize the contribution to our error budget from the FGS Optical Field 
Angle Distortion calibration (Benedict et al. 1999), we held the spacecraft orientation fixed 
for the three prime epochs as shown in Table 1. The last three data sets (epochs 7 — 9 
at a different spacecraft roll) were obtained significantly after the bulk of the observations. 
They maximize the signatures of stellar parallax and proper motion of Gl 876, so that these 
systematic motions could be precisely modeled. Table 1 contains mean epochs for the nine 
distinct observation campaigns. 

We obtained classification spectra (27 October 2000) of our astrometric reference stars 
with the European Southern Observatory Faint Object Spectrograph and Camera (EFOSC) 
on the 3.6m telescope. This effort yields a more accurate absolute parallax for Gl 876, 
reducing the contribution of distance uncertainty to the companion mass determination. 
Reference stars were classified by a combination of template matching and line ratios. The 
half-fringe produced by the FGS walk down to fine-lock (Nelan & Makidon 2002) indicated 
that ref-4 is a binary star with AV ~ 1. With a system brightness, V = 11.02, the primary 
has an intrinsic V = 11.5. Absolute magnitudes for these spectral types from Cox (2000), 
V magnitudes from FGS photometry, and an adopted upper limit Ay = 0.1 from Schlegel 
et al. (1998) at the galactic latitude of Gl 876, produced the absolute parallaxes listed in 
Table 2. These enter our modeling as observations with error. 

With the positions measured by FGS 3 we determine the scale, rotation, and offset 
plate constants for each observation set relative to an arbitrarily adopted constraint epoch. 
The Gl 876 reference frame contains 5 stars. We employ a six parameter model and neglect 
corrections for lateral color discussed in Benedict et al. (1999), because the orientation was 
held constant for the epochs that were expected to yield the largest perturbation signals. 

As in all our previous astrometric analyses, we employ GaussFit (Jefferys, Fitzpatrick, 
& McArthur 1987) to minimize % 2 . The solved equations of condition for the Gl 876 field 
are: 



A, and B are scale and rotation plate constants, C and F are offsets; R x and R y are 
radial terms; fi x and /j y are proper motions; At is the epoch difference from the mean epoch; 
P a and Ps are parallax factors; and ir x and n y are the parallaxes in RA and Dec. We obtain 
the parallax factors from a JPL Earth orbit predictor (Standish 1990), upgraded to version 



f = Ax + By + C + R x (x 2 + y 2 ) - fi x At - P a iv x 
rj = -Bx + Ay + F + R y (x 2 + y 2 ) - fi y At - P s n y 



(1) 
(2) 
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DE405. Orientation to the sky is obtained from ground-based astrometry (USNO-A2.0 
catalog, Monet 1998) with uncertainties of ±0?05. 

We treat reference star proper motions (from the UCAC catalog, Zacharias et al. 2000) 
as observations with error. Results are shown in Table 3. The standard errors resulting 
from the solutions for relative position, parallax, and the final adjusted proper motions 
of Gl 876 and the reference stars are ~ 1 milliarcsecond (mas). The resulting reference 
frame 'catalog' in £ and r] standard coordinates transformed to RA and DEC relative to 
reference star ref-4 was determined with average uncertainties < a^ >= 0.7 and < a v >= 0.5 
mas. Figure 1 shows histograms of target and reference star residuals obtained from our 
astrometric modeling, prior to including the perturbation in the model. 

3. The Astrometric Detection of the Gl 876 Perturbation 

Once proper motion and parallax have been determined, we search for any remaining 
systematic patterns in the Gl 876 residuals. After transforming all residuals to an RA- 
Dec orientation, we plot them (Figure 2) phased to a period of Gl 876b as determined by 
RV measurements (Marcy et al. 2001), P = 61.02 days. Having phased the bulk of the 
astrometric observations to occur near peri- and apastron, the perturbation appears as a 
position difference between these two orbital phases in RA and Dec (Figure 2). The straight 
lines are least-squares fits to all the residuals in both RA and Dec, where the slopes in RA 
and Dec are 5a and la results, respectively. These slopes are determined almost entirely 
by the clumps of residuals at = 0.26 (periastron) and <fi = 0.72. By definition, the total 
vector separation between the residual clumps at peri- and apastron is twice the perturbation 
semi-major axis, 2a. Hence, from astrometry alone we obtain a = 0.3 ±0.1 mas. As in any 
traditional binary star measurement the photocentric orbit as measured by FGS 3 can be 
influenced by the brightness of the secondary (Benedict et al. 2001). Our simulations and 
past measurements of actual binary stars indicate that Am = 5 is the largest Am that 
might shift the primary photocenter by an amount equal to the astrometric precision of our 
final semi-major axis result discussed in Section 4. Assuming this Am and any of the lower 
main sequence mass-luminosity relations referenced below, the system distance, period, and 
primary mass would predict a 54 mas perturbation. As evidenced by the small perturbation, 
the companion must have Am » 5. The barycentric and photocentric position of Gl 876 
are identical within our observational precision. 
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4. The Gl 876b Perturbation Size, Inclination, and Mass from Combined 

Astrometry and Radial Velocities 

Considering only astrometry at peri- and apastron (the two densest clumps seen in 
Figure 2), we could not distinguish between the inclination and orientation of the pertur- 
bation. Properly, our modeling considers the astrometry at all orbital phases. Our analysis 
techniques for defining orbits by combining HST astrometry and RV have been described 
by Benedict et. al. (2001). There we determined all reference and target star parameters 
(including the traditional orbital elements) in a simultaneous solution incorporating both 
astrometry and RV. For Gl 876 we also employ techniques to combine RV from different 
sources (Hatzes et al. 2000) and, additionally, apply several constraints. First, we constrain 
all plate constants to those determined from our astrometry-only solution (equations 1 and 2, 
Section 2). Second, we constrain the orbital elements K, e, P, and uj to those determined by 
us (simultaneously, for both the inner and outer companion) from all available RV data. The 
final constraint is the same one applied in Benedict et al. (2001), the relationship between 
astrometric and RV determined parameters discussed in Pourbaix & Jorissen (2000). 

a sin i PK\sqrt{l — e 2 ) 
itabs ~ (2tt) x 4.7405 

where quantities derived only from our astrometry (parallax, 7r a 5 S , primary perturbation size, 
a, and inclination, i) are on the left, and quantities derivable for the outermost companion 
from RV (the period, P, eccentricity, e, and the RV amplitude for the primary, Ki), are on 
the right. 

With RV measurements complete through 2001 (Keck, Lick, CORALIE, ELODIE), all 
astrometric measurements, and the constraints as described, we solve for parallax, proper 
motion, and the semi-major axis, orbit orientation, and orbit inclination for the perturbation 
caused by the outer companion. The model is aware of the existence of the inner companion, 
and the RV are treated as a sum of the velocities due to both Gl 876b and Gl 876c. But, 
for now the astrometric signature is treated as a perturbation due only to Gl 876b. For the 
parameters critical in determining the mass of Gl 876b we find a parallax, 7r a b s = 214.6 ±0.2 
mas, a perturbation size, a = 0.25 ±0.06 mas, and an inclination, i = 84° ±6°. These, proper 
motion, and the other constrained orbital elements are listed in Tables 4 and 5, along with 
our mass for Gl 876b. The planetary mass critically depends on the mass of the primary 
star, for which we have adopted M* = 0.32 ±O.O5M (Henry & McCarthy 1993, Delfosse et 
al. 2000, Henry & Torres 2002). The first value in Table 5, M b = 1.89 ± 0.34M Jnp , assumes 
no uncertainty in the mass of the primary star. The second value, Mf, = 1.9 ± 0.5Mj up , 
incorporates the present uncertainty in M*. Table 4 also shows satisfactory agreement 
for the parallax and proper motion of Gl 876 from HST, HIPPARCOS (HIP 113020), and 
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the Yale Parallax Catalog (YPC95). If co-planar with Gl 876b, Gl 876c will produce a 
perturbation of order 0.05 mas. 

To confirm the validity of our combined RV and astrometric modeling, we form two 
normal points with the residuals from the data-rich epochs 1, 3, and 4 (Table 1). We can 
collapse the clouds of residuals seen in Figure 2 to normal points, because the astrometric 
residuals are distributed normally (Figure 1). These two best-determined residuals at peri- 
and apastron are plotted in Figure 3. The errors are standard deviations of the means on 
each axis. Comparing with four of the many possible orbits permitted by the equation 3 
constraint, we find satisfactory consistency with the innermost orbit produced by our derived 
parameters (Table 5). 

We have yet to modify our model to treat the mutual interactions of the two components, 
which are expected to result in time- varying orbital elements (Laughlin & Chambers 2001, 
Rivera & Lissauer 2001, Ji, Li, & Liu 2002, Nauenberg 2002). We expect the effect of these 
changes to be small, because the perturbation magnitude is established using observations 
spanning only a single orbital period. As a check we input to our model the dynamically 
determined RV values of Laughlin & Chambers (2001) and obtained a companion mass equal 
to that found by our dual Keplerian approach. We will eventually include evolution of orbital 
elements, which may further reduce the mass uncertainty for Gl 876b. 

GFB, BEM, and WHJ received support through HST NASA grants GO-08102 and GO- 
09233 (T. Forveille, P.I.). We thank Otto Franz, Tom Harrison, and Denise Taylor for useful 
discussions, and an anonymous referee for clarifying suggestions. 
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Table 1. Gl 876 Log of Observations 



Epoch 


JD-2400000 


</> a 


N ofes b 


HST Roll c 


1 


51333.6447 


0.259 


30 


295.87 


2 


51345.1173 


0.447 


5 


295.87 


3 


51361.6886 


0.719 


45 


295.87 


4 


51394.6221 


0.259 


30 


295.87 


5 


51405.9592 


0.444 


5 


295.87 


6 


51696.292 


0.202 


5 


295.87 


7 


51851.7516 


0.75 


5 


98.005 


8 


52229.8838 


0.947 


4 


114.281 


9 


52241.8424 


0.143 


4 


112.886 



a Arbitrary Gl 876b orbital phase, assuming P = 
61.02 days. 

b Number of Gl 876 observations at this epoch 
C FGS 3 orientation 
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Table 2. Astrometric Reference Stars: Spectral Classifications and Spectrophotometric 

Parallaxes 



ID 


SpT 


ya 






7r abs (mas) 


ref-2 


G5 V 


14.47 


5.1±0.4 


0.1 


1.4±0.3 


ref-3 


G8 V 


13.90 


5.5±0.4 


0.1 


2.2±0.4 


ref-4 


GO V 


11.50 d 


4.4±0.4 


0.1 


3.8±0.7 


ref-5 


K2 V 


13.41 


6.4±0.4 


0.1 


4.2±0.8 


ref-6 


Kl V 


15.12 


6.2±0.4 


0.1 


1.8±0.3 



a From FGS photometry reduced and calibrated as 
per Benedict et al. (1998) 

b From Cox (2000) 

c From Schlegel et al. (1998) 

d After removing contribution from a stellar com- 
panion with AV = +1.0 



Table 3. Gl 876 and Reference Star Astrometry 



Star 


RA (arcsec) a,b 


DEC (arcsec) a ' b 


7r a6s (mas) b 


/j, ra (mas) b 


/i DEC (mas) b 


Gl 876 


262.9840±0.0005 


95.3963±0.0003 


214.6±0.6 


954.1±0.8 


- 674.4±0.6 


ref-2 


174.9630±0.0005 


93.9411±0.0005 


1.4±0.1 


5.2±0.9 


- 15.5±0.7 


ref-3 


0.0000±0.0007 


0.0000±0.0008 


2.3±0.2 


13.8±1.2 


3.4±1.1 


ref-4 


-11.0933±0.0008 


-139.6018±0.0005 


3.7±0.2 


-42.3±0.5 


- 43.2±0.4 


ref-5 


663.3341±0.0008 


184.9227±0.0004 


4.2±0.2 


6.3±1.4 


-1.6±1.3 


ref-6 


592.0859±0.0009 


104.4975±0.0006 


1.8±0.1 


-13.9±1.3 


- 3.9±1.2 



a RA and DEC are positions relative to ref-3 (RA = 22 /l 52 m 58 ! :56, Dec = -14°17'24"6, 
J2000.0), epoch = 1999.592 

b After adjustment by the model (equations 1 and 2) 
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Table 4. Astrometry of Gl 876 



Parameter 


Value 


HST study duration 


1.8 y 


number of observation sets 


27 


ref. stars < V > 


13.6 


HST Absolute Parallax a 


214.6±0.2 mas 


HIPPARCOS Absolute Parallax 


212.7 ±2.1 mas 


YPC95 Absolute Parallax 


211.9 ± 5.4 mas 


HST Relative Proper Motion a 


1168.3 ± 1.2 mas y" 1 


in pos. angle 


125?3 ±0?1 


HIPPARCOS Proper Motion 


1174.2 ± 5.4 mas y" 1 


in pos angle 


125? 1 ±0?6 


YPC95 Proper Motion 


1143 mas y _1 


in pos. angle 


123?5 



a Values come from modeling RV and astrometry simul- 
taneously (Section 4) 
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Table 5. Orbital Elements of Perturbation Due to Gl 876b 



Parameter 


Value 


a (mas) 


0.25 ± 0.06 


a (AU) 


0.0012 ± 0.0003 


i 


84°± 6° 


P(days) a 


61.02 ± 0.03 


T (JD) a 


2450107.87 ± 1.9 


e a 


0.10 ± 0.02 


n 


25°± 4° 




338?96 ± 0?36 


Ki (km s _1 


) a 0.210 ± 0.005 


M» (M ) 


0.32±0.05 


M fe (M jMp ) 


1.89±0.34 


M fe (M Jup ) 


1.9±0.5 b 



a Constrained to values deter- 
mined from RV measurements. 
b Error includes M* uncertainty. 
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Fig. 1. — Histograms of x and y astrometric residuals obtained from modeling Gl 876 and 
reference stars with equations 1 and 2. Fits are Gaussians with indicated standard deviations. 



0.0 0.2 0.4 0.6 0.8 1.0 
Arbitrary Phase (P = 61 d .02) 

Fig. 2. — Raw residuals for Gl 876 phased to the Gl 876b companion period, P = 61.02 days, 
determined from radial velocities. The offset between the two primary phases indicates the 
detection of a perturbation. 
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Fig. 3. — Four orbits permitted by the Pourbaix & Jorissen constraint (equation 3) for the 
perturbation of due to Gl 876b. The innermost orbit is determined from the simultaneous 
astrometry-RV solution discussed in Section 4, and is described by the orbital parameters in 
Table 5. Plotted on this orbit are the phases of the astrometric observations, with the two 
primary phases indicated by large circles. Also plotted (+) are the astrometry-only residual 
normal points at 4> = 0.26 (periastron, lower right) and <p = 0.72. These normal points 
are connected to the derived orbit by residual vectors. The insert table shows the variation 
of inclination as a function of perturbation size as mandated by equation 3. Changing 
inclination changes companion mass. 



